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Hypermassive neutron stars (HMNSs) are equihbrium configurations supported against 
collapse by rapid differential rotation and likely form as transient remnants of binary 
neutron star mergers. Though HMNSs are dynamically stable, secular effects such as 
viscosity or magnetic fields tend to bring HMNSs into uniform rotation and thus lead to 
collapse. We simulate the evolution of magnetized HMNSs in axisymmetry using codes 
which solve the Einstein-Maxwell-MHD system of equations. We find that magnetic 
braking and the magnetorotational instability (MRI) both contribute to the eventual 
collapse of HMNSs to rotating black holes surrounded by massive, hot accretion tori and 
coUimated magnetic fields. Such hot tori radiate strongly in neutrinos, and the resulting 
neutrino-antineutrino annihilation could power short-hard GRBs. 
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c/3 ■ Short-hard gamma-ray bursts (SGRBs) emit large amounts of energy in gamma 

rays^ with durations ~ 10~^-2 s and may originate from binary neutron star merg- 
ers. ^'^ If the total mass of the binary is below a certain threshold, a hypermassive 
neutron star (HMNS) likely forms as a transient merger remnant.'''* HMNSs have 



1. Introduction and Methods 



I masses larger than the maximum allowed mass for rigidly rotating neutron stars 

and are supported against collapse mainly by rapid, differential rotation. 

We have performed general relativistic magnetohydrodynamic (GRMHD) simu- 
lations of differentially rotating F-law HMNSs'^''*'^ using two new GRMHD codes. ^'^ 
Here, we consider an HMNS collapse model with a more realistic hybrid EOS. 
For details of the EOS, the chosen HMNS model, and the simulation, see Shibata 
et al.^ and Duez et al.'' We construct a differentially rotating HMNS with mass 
M = 2.65Mq and angular momentum J = 0.82GAP/c. This HMNS is similar to 
one found in a BNS merger simulation.* A small seed poloidal magnetic field is 
introduced into the HMNS as described in Shibata et al.^ These calculations are 
the first in general relativity to self-consistently generate a candidate GRB central 
engine (i.e., a rotating BH surrounded by a magnetized torus) from non-singular 
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Fig. 1. Upper panels: Density contours (solid curves) and velocity vectors at the initial time and 
at a late time. The contours are drawn for p = 10^^ g/cm^ X 10""'*' g/cm^ (i = 0—9). In the 
second panel, a curve with p = 10^^ g/cm"^ is also drawn. The (red) circle in the lower left of the 
second panel denotes an apparent horizon. The lower panels show the poloidal magnetic field lines 
at the same times as the upper panels. 



initial data. 



2. Results and Discussion 

In Figure 1, we sliow snapsliots of the meridional density contours, velocity vec- 
tors, and poloidal magnetic field lines at the initial time and at a late time. The 
differential rotation of the HMNS winds up a toroidal magnetic field, which then 
begins to transport angular momentum from the inner to the outer regions of the 
star (magnetic braking), inducing quasistationary contraction of the HMNS.^ Af- 
ter the toroidal field growth saturates, the evolution is dominated by the MRI,^^ 
which leads to turbulence, thus contributing to the angular momentum transport. 
The star eventually collapses to a BH, while material with high enough specific 
angular momentum remains in an accretion torus. The accretion rate M gradu- 
ally decreases and eventually settles down to M ~ lOAfo/s, giving an accretion 
timescale of Mtorus/-^ ~ 10 ms. 

To explore the properties of the torus, we calculate the surface density E and the 
typical thermal energy per nucleon, u. We find u « 94 MeV/nucleon, or equivalently, 
T « 1.1 X 10^^ K. Because of its high temperature, the torus radiates strongly 
in thermal neutrinos. ^'^'^^ However, the opacity inside the torus is approximately 
K ~ 7 X lQ^^^Ti2 CTP? g^^. The optical depth is then estimated as r kS ~ 
7200I]i8rf2, so that the neutrinos are effectively trapped." Here, Tia = T/IO^^ 
Eis — S/10^^ g cm^^. This regime of accretion has been described as a neutrino- 
dominated accretion flow (NDAF).^^ 

We note that the properties of the torus are not specific to the chosen initial 
data. For example, we consider the model labeled star A in Duez et al.'* This is 
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also an HMNS model which collapses to a BH surrounded by a hot accretion torus 
under the influence of magnetic fields. However, star A has simple F-law EOS (not 
a hybrid EOS) and has a different rotation profile and compactness. In this case, we 

find that the disk has u « 5 MeV/nucIcon, which gives an optical depth of about 
70. Thus, the evolution of star A also produces a hyperaccreting NDAF. 

Returning to the hybrid EOS HMNS model, we estimate the neutrino 
luminosity in the optically-thick diffusion limit. ^® Wc obtain L,^ ^ 2 x 
]^q53 erg/s(i?/10 km)^rf^2^r8^j which is comparable to the neutrino Eddington lu- 
minosity."'^^ A model for the neutrino emission in a similar flow environment with 
comparable gives for the luminosity due to vv annihilation L^p ~ 10^'^ ergs/s.^^ 
Since the lifetime of the torus is ^ 10 ms, the total energy, i^^p ^ 10^* ergs, 
may be suflacient to power SGRBs as long as the emission is somewhat beamed.^'' 
Our numerical results, combined with accretion and jet models, thus suggest 
that magnetized HMNS collapse is a promising candidate for the central engine of 
SGRBs. 
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